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ABSTRACT 

We determine an absolute calibration of the initial mass function (IMF) of early-type galaxies, by 
studying a sample of 56 gravitational lenses identified by the SLACS Survey. Under the assumption 
of standard Navarro, Frenk & White dark matter halos, a combination of lensing, dynamical, and 
stellar population synthesis models is used to disentangle the stellar and dark matter contribution 
for each lens. We define an "IMF mismatch" parameter a =M^ g in /M^^ n as the ratio of stellar mass 

inferred by a joint lensing and dynamical models (M L ) to the current stellar mass inferred from 
stellar populations synthesis models (M^ s in ). We find that a Salpeter IMF provides stellar masses 
in agreement with those inferred by lensing and dynamical models ((log a) = 0.00 ± 0.03 ± 0.02), 
while a Chabrier IMF underestimates them ((log a) = 0.25 ± 0.03 ± 0.02). A tentative trend is found, 
in the sense that a appears to increase with galaxy velocity dispersion. Taken at face value, this 
result would imply a non universal IMF, perhaps dependent on metallicity, age, or abundance ratios 
of the stellar populations. Alternatively, the observed trend may imply non-universal dark matter 
halos with inner density slope increasing with velocity dispersion. While the degeneracy between the 
two interpretations cannot be broken without additional information, the data imply that massive 
early- type galaxies cannot have both a universal IMF and universal dark matter halos. 

Subject headings: gravitational lensing — galaxies: elliptical and lenticular, cD — galaxies: evolution 
— galaxies: formation — galaxies: structure 



1. INTRODUCTION 

The initial mass function (IMF) is a fundamental char- 
acteristic of a simple stellar population. Measuring the 
IMF from resolved stellar populations in the local uni- 
verse hjs_be^ ainajOTastrophy_s^ problem for decades 
(e.g.. ISameteil Il955l : IChabrienl2003h . 

From the point of view of understanding distant unre- 
solved stellar populations, the IMF holds the key to in- 
terpreting observables such as colors and their evolution 
in terms of star formation history and chemical enrich- 
ment history. Among the many astrophysical problems 
where the IMF plays a key role, in this paper we will 
focus on the determination of the stellar mass of galax- 
ies from the comparison of stellar populations synthesis 
models (e.g., iBruzual fc Charlotl [20031: iMarastonl [20051) 
with broad band colors. Stellar masses derived in this 
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manner are often used to construct stellar mass func- 
tions and to st udy the demographics of galaxies over cos- 
mic time ('e.g.jBrinchmann fc Ellis 2000: [Fbntana et alJ 
[20061 IBundv et al.ll2006l )~ 

Unfortunately, the form of the IMF is degenerate with 
the derived stellar masses, since the luminosity is typ- 
ically dominated by stars in a relatively narrow mass 
range. To facilitate comparison between results of dif- 
ferent groups, it is therefore common practice to assume 
a standard, universal, a nd non-evolving IMF. Popular 
choice s are the so-called iSalpeterl ([1955) and iChabrierl 
(2003) IMFs. However, this strategy leaves behind a 
systematic error that is believed to be of the order of 
a factor of two in stellar mass. Furthermore, the er- 
ror need not correspond to a global error on the mean, 
but it could depend on the conditions of star formation 
and therefore for example on galaxy type, environment, 
metallicity and star formation epoch. Recent work, for 
example, has called into question the universality of the 
IMF and suggested it may be evolvi ng with cosmic time 
e.g., van Dokkuml feOOS: Dav1 l2008h . 

The goal of this paper is to determine the absolute nor- 
malization of stellar masses for early-type galaxies, its 
scatter from object to object, and its dependence on sec- 
ondary parameters such as stellar velocity dispersion or 
stellar mass, by combining three independent probes of 
mass. Previous works have attempted to do this by com- 
paring stellar masses determined from stellar populations 
synthesis mo dels with those infer red from gravitational 
lensi ng (e.g.. [Ferreras et al.ll2008|) or stellar kinematic s 
(e.g., iPadmanabhan et al J 12004 iCappellari et al.l [20061) . 
However, the combination of the two latter techniques 
is particularly powerful as it allows one to reduce many 
of the degeneracies and disentangle the stellar and dark 
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component (e.g-. lTreu fe Koopmansl[200l ). 

To meet our goal we exploit the large and homogeneous 
sample of strong gravitat ional lenses discovered by the 
Sloan L enses ACS Survey (IBolton et al.ll2006HTreu et all 
20061: IKoopmans et atl 120061: IGavazzi et alj 120071: 
Bolton et all I2008at IGavazzi et all 120081 : IBolton et al l 
2008bi: IAuger et al.ll2009D . These papers showed that the 
SLACS lenses are statistically indistinguishable within 
the current level of measurement errors from control 
samples in terms of size, luminosity, surface brightness, 
location on the fundamental plane, environment, stellar 
and halo mass. Thus our results can be generalized to 
the overall population o f early -type galaxies. 

Recently. iGrillo et alj (|2009f ) used ground based pho- 
tometry to derive stellar masses for a subset of SLACS 
lenses. They compared the inferred stellar mass fractions 
wi th average tensin g stell ar m ass fractions d etermined 
bv IKoopmans et all (|2006l ) and IGavazzi et all (|2007f ) for 
subsamples of the SLACS lenses, finding a general agree- 
ment. The analysis presented here takes several steps 
forwa rd: i) space-based HST photometry (|Auger et al.l 
2009) extending into the near infrared is used for stellar 
masses; ii) individual stellar fraction estimates from lens- 
ing and dynamical models are used for each galaxy; iii) 
a Bayesian framework is adopted to determine the IMF 
normalization and errors for each galaxies. This progress 
allows us to investigate for the first time the scatter of 
the IMF and its possible dependency on galaxy velocity 
dispersion and hence, e.g., redshift of formation of the 
stellar populations or metallicity. 

We assume a concordance cosmology with matter and 
dark energy density O m = 0.3, Qa = 0.7, and Hubble 
constant Ho==100/ikms _1 Mpc _1 , with h = 0.7 when nec- 
essary. Base-10 logarithms are used. 

2. SAMPLE AND DATA 

The sample analyzed in this paper is composed of 
56 of the 58 early-type lens galaxies identified by 
the SLACS Survey, for which a joint lensing and 
dynamical analysis has been perfo rmed (§ | 2.1|. fol- 
lowing the methods describ ed by iTreu fe Koopmansl 
(|2004) and IKoopmans et alj (120061) . Two of the lenses 
of the IKoopmans et al.l (|2009f ) sample have been ex- 
cluded because the available single-band HST photom- 
etry is not sufficien t to determine reliable stellar masses 
(|Auger et aJ J [2009) . All lenses are classified as definite 
(grade "A") based on the identification of multiple im- 
ages and are successfully modeled as a single mass com- 
ponent. A full description of the SLACS Survey and the 
selection process — together with Hubble Space Tele- 
scope images and measured photometric and spectro- 
scopic paramete rs of all the lense s — is given in the 
SLAC S papers (|Bolton et all 120061 12008at lAuger et all 
[20091) . 

2.1. Dark and luminous mass from lensing and stellar 
kinematics 

The lensing models provides a very accurate and pre- 
cise measurement of the mass contained within the Ein- 
stein Radius (MmJ . As disc ussed at length by, e.g ., 
ITreu fc Koopmansl ([2004D and IKoopmans et all (j2006|) . 
the addition of stellar velocity dispersion information al- 
lows one to disentangle stellar and dark matter, given a 
choice of mass models describing the dark matter halo, 



the stellar component and orbital anisotropy. In practice, 
the joint lensing and dynamical analysis can be decou- 
pled because the uncertainty on the lensing measurement 
is negligible with respect t o that associated to the stel- 
lar velocity dispersion (see lBarnabe et aDl2009l . for joint 
analysis and detailed discussion of the methodology) . We 
thus proceed as follows. First, we determine the total 
mass within the Einstein radius by fitting the lensing 
geometry with a gravitational lens model. The inferred 
lensing mass is determined to a few percent precision, 
independent of the specific form chosen to describe the 
gravitational potential of the deflector (a ssumed to be 
a singular isothermal elli psoid by SLACS; IBolton et al.l 
I2008at I Auger et al] 120091) . We then compute the likeli- 
hood of a family of two-component mass models (stellar 
and dark matter) with respect to the stellar velocity dis- 
persion measured by SDSS and determine the range of 
acceptable solutions by computing the posterior proba- 
bility distribution function. Confidence intervals on each 
individual parameter can be obtained by marginalizing 
over the other parameters. 

For t he purpose of this p a per, w e adopt the so called 
NFW (jNavarro et all fl996l Il99l model for the dark 
matter halo: 

pnM{r n) = . (1) 

(r/r b )(l + (r/r b )r 

In accordance with the CDM picture (e.g. Bullock et 
al. 2001, Maccio et al. 2007) we expect the break radius 
rb to be much larger than the effective and Einstein radii. 
This makes the results insensitive to the precise choice of 
r;,. We fix rb = 30 kpc in agreement with the expected 
value for the average virial mass of the SLACS sample 
(~ 10 13 Mq; Gavazzi et al. 2007). 

Th e lum inous c omponent is des cribed as either a 
[Jaffa (flMl or a iHernquisH ([1990) model, which are 
good simple analytic descriptions of the light profile of 
early- type galaxies, a nd bracket the inner slope of the 
Ide Vaucouleursl (| 19481 ) profile. The orbital anistropy of 
the stars is modeled as a constant (3 

where o~g and ay are the tangential and radial component 
of the pressure tensor. For clarity, we will adopt results 
obtained with Hcrnquist-isotropic models as our default 
(i.e., with [3 = 0). As we will discuss further below, none 
of the results of this paper is changed if Jaffe or moder- 
ately anisotropic models (f3 = ±0.25)- con sistent with in- 
depen dent constraints on anisotropy (e.g. iGerhard et al.l 
l2001h - are considered instead. 

Since the mass within the Einstein radius is fixed by 
the lensing geometry- for a given anisotropy and func- 
tional form of the stellar component- the model has just 
one free parameter: the fraction of stellar mass /* inside 
the cylinder of radius equal to the circularized Einstein 
radius. Thus, for each lens, the lensing and dynamical 
analysis produces the full posterior distribution function 
°f f*i p(/*)j assuming a uniform prior in the interval 
[0,1]. The product of Meid an d /* provides the lens- 
ing+dynamical measurement of the stellar mass inside 
the Einstein Radius M^ g in for each lens galaxy. We note 
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that p(/*), and thus p(M^g in ) are typically fairly asym- 
metric, because of the physical requirement that the stel- 
lar mass fraction be less or equal to unity imposed by the 
prior. For illustration purposes we use the median as our 
best estimator of the quantity, but we use the full dis- 
tribution throughout the analysis. The average median 
value of /* is 80% with a dispersion of 17%, consistent 
with the fact that most of the mass inside the Einstein 
Radius, corresponding on average to half the effective ra- 
dius for the SLACS lenses, is accounted for by the stellar 
component. 

The main goal of this paper is to explore the con- 
straints on the IMF that can be gathered by assuming 
a standard universal NFW profile for the dark matter 
halo. More general forms of the dark matter halo profile 
could also be considered with our formalism, including 
for example the so-called generalized NFW profile, where 
the inner slope is allowed to be a free parameter. The 
gNFW profiled includes steeper profiles as well as pro- 
files with a constant inner core, which are believed to be 
appropriate f or some spiral and low surface brightness 
galaxies (e.g.. lSalucci et alll2007l and r eferences therein! 
as well as some clust ers of galaxies (e.g. lSand et al.|[2b08t 
INewman et al.ll2009t and references therein) . In general, 
allowing the inner slope to be a free parameter results 
in a degenerac y between stellar mass fr action and in- 
ner slope (e.g.. iTreu fe Koopmansl [2004). in the sense 
that steeper inner slopes require less stellar mass to ob- 
tain the same stellar velocity dispersion. The degeneracy 
is best reduced by s patially resolved velocity dispersion 
measurements (e.g.. ITreu fc Koopmansl [2~004j) . However, 
even with a single stellar velocity dispersion measure- 
ment, interesting limits on the inner slope and stellar 
mass fraction can be obtained by marginalizing over the 
other parameter with an appropriate prior. A full anal- 
ysis of the inner slope of dark matter halos is left for 
future work when more accurate and spatially resolved 
velocity dispersion measurements will be available to bet- 
ter constrain the inner slope (e.g., Barnabe et al, 2010, 
in preparation). However, the results presented in this 
paper do not change significantly if a gNFW halo with 
uniform prior on the inner slope is considered instead of 
a simple NFW. 

2.2. Stellar mass from stellar populations synthesis 
models 

The second fundamental ingredient of this work is the 
posterior distribution function for the stellar mass de- 
rived by Auger et al. (2009) by applying stellar popu- 
lations synthesis models to multicolor HST photometry. 
In this paper, we consider as our reference stellar masses 
based on the Bruzual & Chariot (2003) models using an 
informative prior on metallicity taken from the spectro- 
scopic study of Gallazzi et al. 2005 (hereafter the Gal- 
lazzi prior; see Auger et al. 2009 for details). To check 
for possible systematics, we also consider stellar masses 
based on an "ignorant" uniform prior on metallicity (see 
Auger et al. 2009 for detail s), and based on stellar popu- 
lations synthesis models by [Maraston| (|200§). F inally, we 
consid er t wo baseline choi ces of the IMF: the iChabrieil 
(20031) and lSalpeterl (|1955l ) models. 

Taking into account the fraction of light inside the 
cylinder, this method provides an independent measure- 
ment of the stellar mass inside the Einstein radius M^ s in . 



It is important to emphasize that the current stellar mass 
is significantly lower than the stellar mass at zero-age, 
due to mass loss during stellar evolution. For a single 
stellar population with Chabrier IMF, ^50% of the ini- 
tial mass is in the form of gas at 10 billions years of age. 
The fraction is significantly smaller for a Salpeter IMF 
(~30%), or if different prescriptio ns for stellar mas s loss 
are adopted, as discussed, e.g., by iMarastonl (|2005l ). Al- 
though the current stellar mass is the standard quantity 
for this kind of analysis it is clear that our joint lensing 
and dynamics analysis is sensitive to all the mass. There- 
fore, if a fraction of baryons lost during stellar evolution 
were to retain the phase space distribution of their parent 
stars, they would also be counted by the lensing and dy- 
namical analysis towards the component distributed as 
light. Most of the residual gas in elliptical g alaxies is be- 
lieved to be currently in the hot phase, (e.g.. lCiotti et al.l 
1991). The exact phase space distribution of the gas 
lost during stellar evolution depends on the complex in- 
terplay of winds, inflows and outflows, cooling, AGN 
heating, interactions with the environ ment and accretion 
of additional "unprocesse d" gas (e.g.. lPellegrini fc Ciottil 
1998; Pipino et al. 2005). Determining the fate of the 
gas is beyond the scope of this paper. However, X- 
ray observations show that the r esidual gas is a small 
fraction of the stellar mass (e.g. iMathews fc Bri ghcnti 
2003; Humphr ey et al.l [20061 ) and therefore most of the 
gas must be either expelled or recycled into secondary 
episodes of star formation. For simplicity, in this anal- 
ysis, we will consider two extrema that should bracket 
the exact solution. In our default scenario, all gas that 
is not recycled is dispersed and is therefore counted by 
our two component model in the dark matter halo. In 
this case, M^g in needs to be compared with the cur- 
rent M^g S in (including of course stellar remnants such as 
black holes and neutron stars). In the other extreme, all 
gas lost retains the distribution function of the stars and 
is therefore counted by the lensing and dynamical two- 
component model in the stellar component. In this latter 
case, which is effectively an upper limit to the dynamical 
importance of residual gas, M^ g in needs to be compared 
with M^ g in at zero-age. 

3. RESULTS 

A comparison of the two independent determinations 
of stellar mass (M^ ° in and M^ n ) is shown in Fig [TJ 
for four combinations of IMF and lensing and dynamical 
models. The two quantities are tightly correlated, with 
scatter consistent with observational errors. Notice that 
for a Salpeter IMF the points lie on average around the 
identity line, and that changing anisotropy of stellar or- 
bits has very little effects on the inferred M^g in . For a 
Chabrier IMF, in contrast, the current stellar mass un- 
derestimates that inferred from stellar populations syn- 
thesis models. Finally, the zero-age M^ s in is larger than 

M*Ein even for a Chabrier IMF. The trends are robust 
with respect to the choice of stellar population synthe- 
sis models or metallicity priors, as illustrated in Figured 
We note that the data appear to suggest that the relation 
between g in and M^ln is not linear: at low masses 
the data appear to lie above the line indicating the iden- 
tity, while at high masses the data appear to lie below 
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Fig. 1. — Comparison between stellar mass in the cylinder 
of radius equal to the Einstein Radius as inferred from lens- 
ing and dynamical models (x-axis) and that inferred from fit- 
ting stellar populations synthesis models to the observed spectral 
energy distribution (y-axis). The solid line indicates the iden- 
tity. Stellar populations synthesis models by Bruzual & Chariot! 
1 20031) are assumed t ogether with an informative metallicity prior 
I Gallazzi et~alll2005l) . 
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Fig. 2. — As in Figure [T] for different choices of stellar popula- 
tion synthesis models and metallicity priors. Salpeter IMF and 
Hernquist stellar component models are assumed. 

the identity. We will return to this point in Section |3~21 
after we discuss in more detail the overall normalization 
in Section [3~T1 



Hernquist + Salpeter 
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Fig. 3. — Posterior distribution function for the IMF mismatch 
parameter a =M^ E in /M^ s in with respect to a Salpeter IMF. 

3.1. Towards an absolute normalization of the IMF of 
massive galaxies: the "IMF mismatch" parameter 

The goal of this paper is to go beyond a comparison 
and determine the absolute normalization of the IMF for 
each galaxy. This is used to investigate its universality 
in terms of intrinsic scatter and dependency on galaxy 
parameters such as velocity dispersion (hence mass) and 
luminosity. For this purpose, we introduce an "IMF mis- 
match" parameter a =Mj; g in /M^ s in . For each galaxy, 
we determine the posterior distribution function for a 
by combining samples drawn from the posterior distribu- 
tion function for M^g in and from that for M^^ n . The 
resulting posterior distribution samples for a assuming 
Salpeter IMF are shown in Figure [3] for illustration. The 
median values of a for each lens are given in Table [1] 
together w ith other key prope rties of the lens galaxies, 
taken from I Auger et all ()2009f ) and references therein. 

The resulting average values of log a for a variety 
of stellar population synthesis models are summarized 
in Table O The statistical uncertainty on (log a) is 0.03 
dex for any given model. The different choices of dynam- 
ical model (Hernquist vs Jaffe, isotropic vs anisotropic) 
influence the average of log a only at a level of 0.02, 
which can then be neglected in the rest of the discus- 
sion and considered as an additional systematic uncer- 
tainty. Confirming the trends shown in Figure [U the 
statistical analysis shows that a Salpeter IMF tends to 
provide on average a much closer match between M^^j n 

and M^ g in than a Chabrier IMF which appears to pro- 
duce M^; n that are systematically lower than M^g in 
(by 0.25 ±0.03 ±0.02 dex). 

Previous authors used comparisons between indepen- 
dent determinations of stellar masses to sel e ct one IMF 
versus another. For example, iGrillo et al.l (|2009f) con- 
cluded that the Chabrier IMF underestimates stellar 
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masses, and therefore pre ferred a Salpeter I MF, based on 
iBruzual fc Charlotl (|2003f ) models. However. IGrillo et all 
( 2009f ) only used the average dark matter fraction as pub- 
lished for a subset of the SLACS lenses and therefore 
could not make a detailed comparison for each individ- 
ual case. We also n ote that the stellar masses inferred by 
IGrillo et alj (2009) for the galaxies in common with this 
study are systematicall y different than the ones adopted 
here. As discussed by lAuger et al.l (|2009| ). this differ- 
ence is due to a combination of photometry differences, 
and choice of priors for metallicity and age parameters 
(e.g., Grillo et al. 2009 assume constant solar metal- 
licity). However the differences are small and do not 
change the overall normalization, which is in very good 
agreement be tween the two studies. 

In contrast. iCappellari et a l. (2006) used different stel- 
lar populations synthesis models (jVazdekis et al.l fl996f ) 
than the ones adopted here and found that stellar masses 
based on a Salpeter IMF were in some cases too high 
compared to those determined with stellar kinematics, 
reaching the opposite conclusion. However - in presence 
of significant statistical errors - this is actually expected, 
even if the Salpeter IMF were a perfect match to the in- 
trinsic IMF, given that measurements will tend to scatter 
below and above the identity line. Indeed, even with our 
own method, there are systems for which M^j n exceeds 

M*Ein f° r a Salpeter IMF, even though the average of 
the median a is close to unity. Another caveat that must 
be kep t in mind when comparing to the Cappcllar i et al.l 
(2006) study is that their sample extends to significantly 
less massive galaxies than ours (velocity dispersions a 
as low as 60 kins -1 , as opposed to our lower limit of 
approximately 200 kms -1 ). Thus the two samples can 
only be compared directly if the IMF does not depend 
on velocity dispersion or on galaxy mass. 

Regardless of the interpretation in terms of a specific 
IMF - which depends also on t he uncertainties of the 
stellar populations models (e.g. iMarastorJ I2005D - we 
emphasize that our method provides an absolute cali- 
bration of the stellar mass. The current stellar masses 
given in Auger et al. (2009) assuming a Salpeter IMF 
multiplied by a are absolutely calibrated against those 
inferred b y the lensing and d ynamical models. On av- 
erage, the lAuger et al.1 (|2009l ) Salpeter masses are cal- 
ibrated to within 0.04 dex, even without applying the 
IMF mismatch parameter. In general, the a values given 
in Table Q] can be used to calibrate any stellar population 
synthesis model, for any arbitrary choice of IMF, and des- 
tiny of the gas lost during stellar evolution. Remarkably, 
the data are consistent with very little intrinsic scatter in 
log a. The upper limit on the intrinsic scatter is 0.09 dex 
(95% CL), i.e. the absolute normalization of the IMF is 
uniform to better than 25%. 

3.2. Universal or not? Trends with galaxy properties 

Within the class of massive early-type galaxies, the 
SLACS lenses span approximately a factor of two in ve- 
locity disper sion and a factor o f ten in luminosity and 
stellar mass (| Auger et al.l I2009T ) . In turn, these quanti- 
ties correlate with the average epoch of formation of their 
stellar populations, as well as their ave rage metalliticity , 
abundance rat ios, and gas content (e.g.,lTreu et al.ll2005l: 
iThomas et al.l 120051 : iGallazzi et al.1 120051 : Uimenez et al.1 



[20071 IPipino"eTaL1l2009l : [Graves et al.ll2009f ). If the IMF 

were to evolve during the epoch of formation of most of 
the stars of early-type galaxies, or if it were to depend on 
the mode of star formation or on the physical condition 
of the progenitor gas, we would expect a to vary across 
our sample. 

To test for signs of mass dependency of the IMF nor- 
malization, we checked for a correlation between a and 
three indicators of galaxy "mass": i) ctsie, the velocity 
disper s ion of the best fittin g lensing model (jBolton et al.l 
2008a- lAuger et al.1 120091 ): ii) cr», the stellar velocity 
dispersion within the SDSS fiber aperture; iii) the to- 
tal V band luminosity c orrected to a comm on redshift 
z = 0.2 as described by lAuger et~aT1 (|2009h . The first 
two choices are motiva ted by several lines of evidence 
(e.g. IGraves et al.l 120091 ) that indicate that velocity dis- 
persion is the most important parameter in determin- 
ing stellar populations. The first quantity correlates well 
with stellar velo city dispersion an d is measured much 
more accurately (|Treu et al.l 120061) . The errors on ctsie 
are effectively negligible with respect to those on stellar 
velocity dispersion, which dominate the error on M^g in . 
Thus, this choice makes the covariance between a and 
ctsie negligible. The canonical stellar velocity dispersion 
cr» suffers from a larger covariance with a due to its larger 
errors. The third quantity, the V-band luminosity, is an 
inferior galaxy mass proxy - because it is sensitive to rel- 
atively minor recent episodes of star formation - and is 
inversely covariant with a because to first order M;^ n is 
proportional to Ly. 

The results are shown in Figure H] using the Salpeter 
IMF as baseline and for our standard stellar populations 
models (adopting a Chabrier IMF would move all the 
points upwards by ~ 0.25 dex, while all the other choices 
introduce negligible changes). A trend with non-zero 
slope is detected for ctsie and cr*. No significant slope 
is found for Ly- The best fit linear relations are found 
to be, with no evidence of intrinsic scatter: 

loga = (1.20 ±0.25) log ctsie -2. 91 ±0.02, (3) 
loga= (1.31 ±0.16) log ct* -3.14 ±0.01, (4) 

log a = (0.11 ± 0.08) logLy + 0.00 ± 0.02, (5) 

where velocity dispersions are in units of kms -1 and Ly 
is in units of 10 11 L y. 

In summary, log a appears to increase with velocity 
dispersion by an amount comparable to the difference 
between a Ghabrier and Salpeter IMF over the range 
probed. The slope of the correlation with luminosity is 
significantly smaller then expected given the correlations 
with velocity dispersion and the Faber Jackson relation. 
The inverse covariance mentioned above is not sufficient 
to explain the discrepancy unless there is significant in- 
trinsic scatter. This may therefore be another indication 
that velocity dispersion and not luminosity is the main 
para meter controlling stellar populations, inc luding the 
IMF (jBernardi et alj|2007t IGraves et al.ll2009T) . 

4. DISCUSSION 
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Fig. 4. — Template mismatch parameter a =M^ jjin/^f^Ein ^ or Salpeter IMF as a function of lensing velocity dispersion (left), stellar 
velocity dispersion (center) and V-band luminosity corrected to z = 0.2. A tentative positive trend with velocity dispersion is observed 
(solid line). The dashed line represents the trend expected for a universal Salpeter IMF. 



Before attempting to interpret our perhaps surpris- 
ing findings, it is important to emphasize a number of 
caveats: i) the sample is relatively small and with a se- 
lection function that strongly favors high velocity dis- 
persion galaxies; ii) the quantities on the two axis of 
Figure 2] are not independent, even though the known 
errors are small enough for covariance not to be caus- 
ing the observed trends; iii) the stellar masses are typi- 
cally based on three-band photometry and therefore we 
can only probe simple star formation histories. For all 
these reasons the results reported here must be consid- 
ered as tentative until verified by larger samples, span- 
ning a broader range of properties, and with the help of 
spectroscopic diagnostics of stellar populations. 

Keeping these caveats in mind, we now discuss the two 
main results of this paper. The first result is that the 
average absolute normalization of the IMF inferred by 
our study is higher than those commonly assumed when 
deriving masses for distant galaxies using colors. Those 
"lighter" IMFs have been preferr ed on the basis of studies 
of loc al stellar populations (e.g. iKroupal [20011 : iChabrierl 
l20f)l as well as on the basis of dynamical argume nts ap- 
plied to spiral galaxies (e.g. iBell fe de Jong]|200ll ). How- 
ever, those measurements do not necessarily apply to the 
stellar populations of massive early-type galaxies, if the 
IMF is not universal, but depends, for example, on metal- 
licity or other c onditions that vary with cosmic time (e.g. 
lElmegreenll2008L and references therein) . The second re- 
sult is the trend in IMF normalization with galaxy ve- 
locity dispersion. Taken at face value, this trend would 
imply that whereas a "light" IMF such as Chabrier's is 
appropriate for systems with a ~ 200 kms -1 (and there- 
fore more or less consistent with the standard conclu- 
sions for Milky- Way- type galaxies and spirals), for the 
most massive systems there is a higher abundance of low 
mass stars. This appears to be quite a dramatic change 
over a factor of ~ 2 in velocity dispersion, correspond- 
ing approxima tely to 0.15 dex chan ge in [Fe/H] and 0.1 
dex in [a/Fe] (Bcrnard i et al.ll2006f ). Part of this trend 
could also be ascribed to larger fraction of gas loss during 
stellar evolution being retained by higher velocity dis- 
persion systems. However, the amount of retained gas 
would have to be comparable to the stellar mass in the 
central regions in order to explain the trend completely, 



which seems to be ruled ou t by X-ray observations (e.g. 
Mathews & Brighenti 2003) . Deep X-ray observations of 
the SLACS sample would be useful to verify exactly how 
much gas is left. More and better spectroscopic data 
are needed to investigate whether this trend is real, or 
whether there are unknown systematic effects at play. If 
the trend were to be confirmed, it would have far reaching 
implications for the determination of the evolving mass 
function of galaxies, because it may change its shape as 
well as its overall normalization. 

Alternatively - if the IMF normalization were indeed 
universal over the mass range spanned by the SLACS 
sample - our finding would imply that one of our as- 
sumptions in deriving g in is not warranted. As we 
have shown, anisotropy or changes in the assumed stel- 
lar mass density profile are not going to be sufficient, as 
they only change a by a few hundredths of a dex at most. 
Only changing our assumed dark matter density profile 
systematically with a would have a sufficiently large ef- 
fect to affect the observed trend. The assumed break ra- 
dius (and hence concentration parameter) has only min- 
imal effects. The dominant parameter in determining /» 
is the inner slope of the dark matter density profile 7: 
stee per halos require a smaller stellar mass fraction (see , 
e.g.. iTreu fe Koopmansl l2002t iKoopmans fc Treul [20031 : 
iTreu fc Koopmansl |2004[ ) . Thus - in order to keep a 
constant and consistent with a Chabrier normalization - 
dark matter halos would need to be NFW-like at the low 
end of the velocity dispersion range and become steeper 
towards the high end. In contrast - in order to keep a 
constant and consistent with a Salpeter normalization - 
dark matter halos would need to be NFW-like at the high 
end of the velocity dispersion range and become flatter 
towards the low end, tending towards an inner constant 
core. 

What could cause the dark matter inner slope to 
steepen with velocity dispersion? Baryons are the pri- 
mary suspect, since this trend is not observed in dark 
matter only simulations. Baryons are indeed dominant 
over this radial range, and they could be responsible 
for making the dark matter density profile steeper tha n 
NFW (e.g. lQnedin etld1l200i Ijiang fc Kochanekll2007f) . 
This baryonic compression would need to be more ef- 
fective for higher velocity dispersion objects in order to 
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explain the observed trend 8 . By steepening the dark 
matter density profile, baryons would also effectively in- 
crease the dark matter fraction within a fixed aperture. 
This in turn would imply a correlation between dark mat- 
ter fraction and velocity dispersion. 

5. CONCLUSIONS 

We combined three independent diagnostics of mass 
(lensing, dynamics and stellar populations synthesis 
models) to determine an absolute normalization of the 
IMF for a sample of 56 early-type galaxies, spanning over 
a decade in stellar mass and a factor of two in velocity 
dispersion, under the assumption of a NFW dark matter 
density profile. On average, the absolute IMF normal- 
ization is found to be close to that of a Salpeter IMF 
and larger than that for a Chabrier IMF. Using the pre- 
scription outlined in this paper, stellar masses based on 
any stellar population synthesis models can be absolutely 
calibrated to better than 20%, a significant progress with 
respect to the range of a factor of ^2 spanned by stan- 
dard choices of the IMF. 

A tentative trend of IMF mismatch parameter 
a=M^g in /Mf^ n with galaxy velocity dispersion is 
found. Two possible explanations, not necessarily mu- 
tually exclusive, are suggested for the observed trend: 

• The IMF is not universal, but rather depends on 
parameters such as metallicity, age, and abun- 
dance ratios of the stellar populations. In order 
to fully explain the observed trend, the normaliza- 
tion of the IMF and thus the abundance of low 
mass stars, must increase from Chabrier-like for 
fj ~ 200 kins" 1 to Salpeter-like for the most mas- 
sive early-type galaxies. 

• Dark matter halos are not universal. For a uni- 
form Chabrier IMF, the observed trend of a with 
a could be explained if the inner slope of the dark 
matter halo were systematically steeper than NFW 
for the high velocity dispersion systems. For a 
uniform Salpeter IMF, the the dark matter halos 
would have to be NFW at the high mass end and 
flatter at lower masses. 

In conclusion, the data are inconsistent with both a 
universal IMF and universal NFW dark matter halos over 
the mass range probed by the SLACS sample. There is a 
fundamental degeneracy between the two interpretations 
that cannot be broken with the current dataset. 

Further tests and more work are required to verify and 
extend our perhaps surprising results. Firstly, we need 

8 While t his manuscrip t was under revision, we became aware of 
the work by Schulz ct al. (2009) who find a similar trend based on 
a joint weak lensing and dynamical of early- type galaxies selected 



to extend our samples to cover a wider range of redshifts, 
masses, and morphological types and thus probe a larger 
variety of IMFs. Secondly, we need to use spectral stellar 
population diagnostics to obtain independent constraints 
on the stellar mass to light ratios as well as on the phys- 
ical parameters that may correlate with IMF normaliza- 
tion. Thirdly, we need to improve the constraints on 
the inner slope of the dark matter halo as a function 
of velocity dispersion to break the current degeneracy. 
At the level of individual galaxies, some progress can be 
achieved by constructing more sophisticated l ensing and 
dynamical models fe.g.. [Barnabe et al.l 12009). In fact, 
the models presented here only use a single measure- 
ment of stellar velocity dispersion from SDSS and the 
total mass enclosed by the Einstein Radius, while more 
radial information can be extracted from both diagnos- 
tics. At the level of joint analysis of sub samples of galax- 
ies, we need to have enough objects so that they can be 
binned by velocity dispersions to perform a weak lensing 
analysis. The addition of weak lensing to the strong lens- 
ing, dynamics, and stellar populations diagnostics would 
allow us to probe systematic variations with velocity dis- 
persion of the dark matter halo shape and of the stellar 
to virial mass to light ratio . If the amount of contraction 
is an important ingredient of the observed trend with ve- 
locity dispersion, we expect to see a parallel trend in the 
overall efficiency of converting baryons into stars, or per- 
haps in the spatial concentration of the stellar component 
relative to that of the halo. 
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TABLE 1 

Basic parameters of the lens galaxies 



ID 




<r* asm L v (z=0.2) 


log a 






(kms" 1 ) (kms" 1 ) (10 11 M ) 





J0029-0055 

J0037-0942 

J0044+0113 

J0216-0813 

J0252+0039 

J0330-0020 

J0728+3835 

J0737+3216 

J0822+2652 

J0841+3824 

J0912+0029 

J0935-0003 

J0936+0913 

J0946+1006 

J0955+0101 

J0956+5100 

J0959+4416 

J0959+0410 

J1016+3859 

J1020+1122 

J1023+4230 

J1029+0420 

J1032+5322 

J1103+5322 

J1106+5228 

J1112+0826 

J1134+6027 

J1142+1001 

J1143-0144 

J1153+4612 

J1204+0358 

J1205+4910 

J1213+6708 

J1218+0830 

J1250+0523 

J1402+6321 

J1403+0006 

J1416+5136 

J1420+6019 

J1430+4105 

J1432+6317 

J1436-0000 

J1443+0304 

J1451-0239 

J1525+3327 

J1531-0105 

J1538+5817 

J1621+3931 

J1627-0053 

J1630+4520 

J1636+4707 

J2238-0754 

J2300+0022 

J2303+1422 

J2321-0939 

J2341+0000 



0.2270 
0.1955 
0.1196 
0.3317 
0.2803 
0.3507 
0.2058 
0.3223 
0.2414 
0.1159 
0.1642 
0.3475 
0.1897 
0.2219 
0.1109 
0.2405 
0.2369 
0.1260 
0.1679 
0.2822 
0.1912 
0.1045 
0.1334 
0.1582 
0.0955 
0.2730 
0.1528 
0.2218 
0.1060 
0.1797 
0.1644 
0.2150 
0.1229 
0.1350 
0.2318 
0.2046 
0.1888 
0.2987 
0.0629 
0.2850 
0.1230 
0.2852 
0.1338 
0.1254 
0.3583 
0.1596 
0.1428 
0.2449 
0.2076 
0.2479 
0.2282 
0.1371 
0.2285 
0.1553 
0.0819 
0.1860 



229±18 
279±14 
266±13 
333±23 
164±12 
212±21 
214±11 
338±17 
259±15 
225±11 
326±16 
396±35 
243±12 
263±21 
192±13 
334±17 
244±19 
197±13 
247±13 
282±18 
242±15 
210±11 
296±15 
196±12 
262±13 
320±20 
239±12 
221±22 
269±13 
226±15 
267±17 
281±14 
292±15 
219±11 
252±14 
267±17 
213±17 
240±25 
205±10 
322±32 
199±10 
224±17 
209±11 
223±14 
264±26 
279±14 
189±12 
236±20 
290±15 
276±16 
231±15 
198±11 
279±17 
255±16 
249±12 
207±13 



217.3±5.0 
285.2±6.6 
268.7±6.2 
347.6±8.0 
234.7±5.4 
251.7±5.8 
256.4±5.9 
292.5±6.7 
270.8±6.2 
247.9±5.7 
346.2±8.0 
360.8±8.3 
242.8±5.6 
283.5±6.5 
223.8±5.1 
318.0±7.3 
253.6±5.8 
215.8±5.0 
253.2±5.8 
303.7±7.0 
267.1±6.1 
208.6±4.8 
249.6±5.7 
217.4±5.0 
239.2±5.5 
314.4±7.2 
242.4±5.6 
254.3±5.8 
285.5±6.6 
220.0±5.1 
253.9±5.8 
285.2±6.6 
251.4±5.8 
254.0±5.8 
243.5±5.6 
293.7±6.8 
224.8±5.2 
287.0±6.6 
204.0±4.7 
336.9±7.7 
235.8±5.4 
256.1±5.9 
207.0±4.8 
221.9±5.1 
317.9±7.3 
281.4±6.5 
222.3±5.1 
284.6±6.5 
273.8±6.3 
311.1±7.2 
247.2±5.7 
238.4±5.5 
300.8±6.9 
289.7±6.7 
259.1±6.0 
262.0±6.0 



0.69±0.05 
1.09±0.06 
0.62±0.06 
1.81±0.12 
0.49±0.04 
0.74±0.07 
0.82±0.05 
1.51±0.09 
0.87±0.06 
0.96±0.09 
1.47±0.07 
2.06±0.20 
0.83±0.06 
0.66±0.04 
0.17±0.01 
1.12±0.08 
0.86±0.06 
0.23±0.01 
0.51±0.04 
0.94±0.07 
0.63±0.04 
0.33±0.03 
0.22±0.02 
0.56±0.04 
0.47±0.03 
0.90±0.06 
0.52±0.04 
0.68±0.04 
0.94±0.08 
0.39±0.03 
0.44±0.02 
0.96±0.05 
0.68±0.06 
0.78±0.07 
1.08±0.06 
1.05±0.05 
0.50±0.03 
0.73±0.05 
0.31±0.02 
1.28±0.10 
1.10±0.09 
1.00±0.07 
0.24±0.01 
0.57±0.03 
1.72±0.16 
1.02±0.08 
0.36±0.03 
1.03±0.06 
0.79±0.05 
1.04±0.05 
0.73±0.05 
0.56±0.03 
0.72±0.04 
1.03±0.05 
0.78±0.05 
0.89±0.06 



-0.07 
0.09 
0.08 
0.09 
-0.50 
-0.06 
-0.18 
0.07 
0.04 
-0.01 
0.13 
0.19 
-0.08 
0.20 
0.03 
0.16 
-0.10 
0.04 
0.01 
0.03 
0.06 
-0.13 
0.11 
-0.17 
-0.06 
0.17 
0.01 
-0.04 
0.23 
-0.01 
0.08 
0.14 
0.17 
-0.01 
-0.13 
0.05 
-0.11 
0.01 
-0.13 
0.11 
-0.41 
-0.01 
-0.03 
0.02 
-0.11 
0.11 
-0.04 
-0.02 
0.06 
0.01 
-0.05 
-0.07 
0.15 
0.08 
0.01 
-0.24 



+0.14 
0.15 
+0.07 
0.07 
+0.11 
0.11 
+0.10 
0.15 
+0.26 
0.38 
+0.15 
0.20 
+0.14 
0.16 
+0.08 
0.09 
+0.13 
0.13 
+0.17 
0.16 
+0.08 
0.09 
+ 0.20 
0.44 
+0.12 
0.12 
+0.14 
0.16 
+0.10 
0.12 
+ 0.08 
0.09 
+0.14 
0.17 
+0.08 
0.10 
+0.12 
0.12 
+0.13 
0.13 
+0.13 
0.14 
+0.11 
0.11 
+ 0.11 
-0.11 
+0.10 
0.12 
+0.06 
0.06 
+0.09 
0.09 
+0.11 
0.12 
+0.15 
0.27 
+ 0.11 
0.10 
+0.13 
0.12 
+0.07 
0.07 
+0.07 
0.08 
+0.10 
0.08 
+0.11 
0.13 
+0.08 
0.07 
+0.09 
-0.12 
+ 0.11 
-0.13 
+0.15 
0.20 
+0.06 
0.07 
+0.14 
0.18 
+0.31 
0.48 
+0.14 
0.21 
+0.06 
0.07 
+0.08 
-0.09 
+0.21 
-0.37 
+0.10 
-0.08 
+ 0.11 
0.14 
+0.15 
0.22 
+0.09 
0.09 
+ 0.10 
0.11 
+0.10 
0.11 
+0.12 
0.14 
+0.08 
0.10 
+0.10 
0.13 
+0.08 
0.09 
+0.19 
0.26 



Note. — The IMF mismatch parameter a is given with respect to a 
Salpeter IMF assuming BC03 models and Gallazzi metallicity prior. cr» 
is the SDSS-measured stellar velocity dispersion within th e spectroscopic 
aperture, as given in SLACS paper V iBolton ct al. 20083). 
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TABLE 2 
Average IMF mismatch parameter 
((log a)) 



IMF SSP models Z Prior (log a) 



Salpctcr BC03 

Chabricr BC03 

Salpctcr BC03 

Chabricr BC03 

Salpctcr M05 

Salpctcr M05 



Gallazzi 0.00 

Gallazzi 0.25 

Uniform 0.03 

Uniform 0.27 

Gallazzi 0.05 

Uniform 0.06 



Note. — 



Statistical errors arc 0.03 dcx. 



